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ABSTRACT 



Context. Significant quantities of magnetized plasma are transported from the Sun to the interstellar medium via Interplanetary 
Coronal Mass Ejections (ICMEs). Magnetic Clouds (MCs) are a particular subset of ICMEs, forming large scale magnetic flux 
ropes. Their evolution in the solar wind is complex and mainly determined by their own magnetic forces and the interaction with the 
surrounding solar wind. 

Aims. In this work we analyze the evolution of a particular MC (observed on March 1998) using in situ observations made by two 
spacecraft approximately aligned with the Sun, the first one at 1 AU from the Sun and the second one at 5.4 AU. We study the MC 
expansion, its consequent decrease of magnetic field intensity and mass density, and the possible evolution of the so-called global 
ideal-MHD invariants. 

Mefhods. We describe the magnetic configuration of the MC at both spacecraft using different models and compute relevant global 
quantities (magnetic fluxes, helicity and energy) at both helio-distances. We also track back this structure to the Sun, in order to find 
out its solar source. 

Results. We find that the flux rope is significantly distorted at 5.4 AU. However, we are able to analyze the data before the flux rope 
center is over-passed and compare it with observations at 1 AU. From the observed decay of magnetic field and mass density, we 
quantify how anisotropic is the expansion, and the consequent deformation of the flux rope in favor of a cross section with an aspect 
ratio at 5.4 AU of « 1.6 (larger in the direction perpendicular to the radial direction from the Sun). We quantify the ideal-MHD 
invariants and magnetic energy at both locations, and find that invariants are almost conserved, while the magnetic energy decays as 
expected with the expansion rate found. 

Conclusions. The use of MHD invariants to link structures at the Sun and the interplanetary medium is supported by the results of 
this multispacecraft study. We also conclude that the local dimensionless expansion rate, that is computed from the velocity profile 
observed by a single spacecraft, is very accurate for predicting the evolution of flux ropes in the solar wind. 

Key words. Magnetic fields, Magnetohydrodynamics (MHD), Sun: coronal mass ejections (CMEs), Sun: magnetic fields. 
Interplanetary medium 



1. Introduction 

Coronal Mass Ejections (CMEs) are explosive events that re- 
lease energy in the solar atmosphere. The interplanetary coun- 
terparts of CMEs are solar wind (SW) structures known as 
interplanetary coronal mass ejections (ICMEs). Among them 
there is a subset, called magnetic clouds (MCs), which exhibit a 
smooth rotation of the magnetic field direction through a large 
angle, enhanced magnetic field strength, low proton tempera- 
ture and a low proton plasma beta, /3p. MCs are formed by 
large scale magnetic flux ropes carrying a large amount of mag- 
netic helicity, magnetic flux and energy away from the Sun. The 
main characterist i cs of t hese structures have been enumerated by 
iBurlaga & Kleinl ([T980h . 

Several authors have consider MCs as static flux ropes 
(see, e.g ., Goldstein 1983; Burlaga 1988; Lepping et al] |1990l: 
iBurlagal 1 19951: iLvnch et al.l l2003l) . Their ' magnetic fields 
have been freque ntly modeled using the Lundquist's model 
dLundquisj Il950l) . which considers a static and axially- 
symmetric linear force-free magnetic configuration. Many de- 
viations from this mod el have been also st udied: e.g., non- 
linear force-free field s ( FarrugiaetalJ 19991), non - force free 
fields (iMuUigan et alJ Il999bt iHidalgo et al.1 l200l ICid et al.1 



2002) 



2001 



and several non-cylindric al models ( Hu & Sonnerud 
IVandas & Romashetsl 120021: iDemouUn & Dassoll2009br 



all of them being static. These models are recurrently used to fit 
in situ magnetic field measurements within MCs to reconstruct 
the whole flux rope structure. These techniques have also been 
tested by replacing the observations by the local values found 
in a numerical simulation, and the output of the models has 
been compared to the known original full simulation (iRilev et al.l 
|2004|). The results of these comparisons show that these in situ 
techniques can reproduce relatively well the magnetic structures 
when the spacecraft is crossing the MC near its main axis. 

In many cases, MCs present clear characteristics of ex- 
pansion (e.g. iLepping et al.l I2003L l2008h . so several dynami- 
cal models have been developed to describe these clouds dur- 
ing their observation time. Some of these flux rope models 



dParrugia et al. 19931 Osherovich et a 


IJ 1993bl: Farrugia et all 


11997: IShimazu & Marubashil 20001 


Nakwacki et alj l2008bl). 



while other models include expansion in both directions, radial 
and axial dShimazu & Vandas 2002; Berdichevs ky et al. 20031; 
iDemouUn & Dassol2o69al; iNakwacki'et al.l2008al). A dynamica l 
model with an elliptical shape was derived bv lHidalgd (l2003h . 
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while a model of the expansion with an anisotropic self-similar 
expansion in three orthogonal directions was worked out by 
[Demoulin et al. (2008). 

From single spacecraft observations we cannot directly infer 
the global structure of the flux ropes and their evolution through 
the interplanetary medium because they are one single point lo- 
cal measurements. Several strategies are used to derive more in- 
formation on MCs with multi-spacecraft data, as follows. 

With two spacecraft located at a similar distance from the 
Sun and separated by a distance of the order of the cross-section 
of the encountered flux rope, the in situ observations provide 
data at different parts of the flux rope being negligibly affected 
by the evolution. This is used to test the technique comput- 
ing the magnetic field in the cross section from the data of 
one spacecraft and/or to have a more accurate reconstruction 
of flie magnetic fie ld (Mulligan & Russell 2001; Liu et al. 200l 
iKilDua et al.ll2()09l:lMostret alllMOQl) . 

When the two spacecraft positions are viewed from the Sun 
with a significant angle, typically in the interval [10°, 80°], one 
can usually derive an estimation of the extension of the flux 
rope, or at least an estimation of the extension of the pertur- 
bation (e.g. the front shock) induced by t he propagation of 
the flux rope in the interplanetary medium (|Cane et al]|1997l : 



iMulUgan & Russelll2001l:lReisenfeld et al.l2003h . A larger num- 
ber of spacecraft permits to constrain more the evolving mag- 
netic s tructure, such as in the case analyzed by iBurlaga et al.l 
([1981 . Such studies have been extended to cases where the 
spacecraft are well separated in solar distance with on e space- 
craft near Earth and the other one at a few AUs ([H amm ond et all 
l^; 'Gosling et al. T99l iLiu et al.ll2006t iFoujlon et al J 120071: 



Rodriguez et al. 2008.) . They show that large MCs/ICMEs have 
large scale effects in the heliosphere (e.g. both at low and high 
latitudes). 

When the two spacecraft are separated by spatial scales of 
the order of one or several AUs, the above analysis should take 
into account the evolution of the MC with solar distance. This 
also implies a more difficult association of the in situ observa- 
tions at both spacecraft. Numerical simulations are then use- 
ful tools to check if the eve nts observed at ea ch spacecraft are 
in fact a unique event (e.g., iRilev et al1l2003h . Since MCs are 
moving mostly radially away from the Sun, the radial alignment 
(line-up) of two spacecraft is a major opportunity to study the 
radial evolution of a MC, as the MC is crossed at a similar lo- 
cation by the two spacecraft. However, it is not common to find 
events observed by two nearly radially aligned spacecraft. One 
case was observed by Helios- 1,2 close to 1 AU and later by 
Voyager- 1,2 at 2 AU, with an angular separation from the Sun 
of about 10° (rest ). 23°) between HeUos-1 (resp. Helios-2) an d 
both Voyager- 1,2 dBurlaga et al.lll981l;IOsherovich et al.]|1993ah . 
A second case was observed by Wind and NEAR spacecraft with 
an angular separation from the Sun of abou t 1° and a ratio of 
solar distances of 1.2 ( iMulligan et al.|[T999al) . A third case was 
observed by ACE and NEAR spacecraft with an angular separa- 
tion from the Sun of a bout 2° and a ratio of solar distances of 1 .8 
(iMulligan et al.ll200Tl) . Finally, a fourth case was observed by 
ACE and Ulysses spacecraft with an angular separat ion from the 
Sun of about 6° and a ratio of solar distances of 5.4 (Sk oug et al.l 
l2000l:lDu et al.l2007l) . This last case has the advantage of a larger 
radial separation, so that the evolution has a larger effect. This is 
the MC selected for a deeper study in this paper. 

From the in situ data at different solar distances of the same 
MC one can infer directly the evolution of the magnetic field and 
plasma quantities. Such radial evolution is otherwise available 
only from a statistic analysis of a large number of MCs observed 



individually at various di s tances (Liuet al. 20051: IWang et aH 
I2005at iLeitner et al l l2007t IGulisano et al.ll201(£ " with possible 
bias coming from the selection of MCs with different properties. 
Another application of line-up spacecraft is to derive the evo- 
lution of global magnetohydrodynamic (MHD) quari tities, such 
as magnetic flux and magnetic helicity (lDassoll2009l) . They are 
main quantities to test if the flux rope simply expands or if a sig- 
nificant part reconnects with the SW field. These global quan- 
tities permit also a quantitative link to the related solar even t 
(iMandrini et al.ll2005[ iLuoni et al.ll200l iRodriguez et'aDl2008l) . 
and set constrain ts on the physical mechanism of the associ- 
ated CME launch (IWebb et al.l2000l:lAttrill et al.l2006l:loiu"etal1 
12007 ft . 

In this paper we further analyze the evolution of a MC 
observed at two different helio -distances, 1 and 5.4 AU 
(ISkoug et al.ll2000l: iDu et al.ll2007l) . This MC has been selected 
because this is to our knowledge the best line-up observations 
of a MC between ACE and Ulysses spacecraft. The observations 
are summarized in Section|2] The velocity and magnetic models 
used to complement the observations are described in Section|3] 
This spacecraft line-up is an opportunity to follow the evolution 
of the flux rope and, in particular, the global MHD quantities 
such as magnetic helicity and flux (Section |4]i. We relate this 
MC to its solar source in Section |5] This complements our un- 
derstanding of the magnetic field evolution. Finally, we discuss 
our results and conclude in Section|6] 



2. Observations 

2.1. Instruments and spacecraft 

We analyze data sets for SW plasma and magnetic field 
from ACE and Uly sses spacecraft. W e use the Magnetic Field 
Experiment (MAG, ISmith et^ll998h with a temporal cadence 
of 16 seconds and the Solar Wind Ele ctron Proton Alpha 
Monitor (SWEPAM. lMcComas et al.lll998l) with a temporal ca- 
dence of 64 seconds for ACE spacecraft. For U lysses space- 
craft, we use Vector Helium Magnetometer (VHM. lBalogh et aH 
11992) for magnetic field observations with a temporal cadence 
of 1 second and Solar Wind Obse rvations Over the Poles of the 
Sun (SWOOPS. lBame et alJll992l) for plasma observations with 
a temporal cadence of 4 minutes. 

When the MC passed through Earth (March 5, 1998) ACE 
was located at a: 1 AU in the ecliptic plane, and in a longitude 
of 164° in the Solar Ecliptic (SE) coordinate system. When the 
cloud was observed by Ulysses (March 25, 1998), this spacecraft 
was located at 5.4 AU from the Sun and very near the ecliptic 
plane, in particular it was at a latitude of 2° and at a longitude of 
158°, in the SE system. Thus, the position of both spacecraft dif- 
fers in 2° for la titude and in 6° for longitude (ISkoug et alJl2000l: 
IDu et al.ll2007l) . This angular separation corresponds to a separa- 
tion distance (perpendicular to the radial direction to the Sun) of 
a 0.6 AU at the location of Ulysses. This very good alignment 
between the Sun and both points of observation of the same ob- 
ject, gives us a unique opportunity to observe the same MC at 
two different evolution stages in the heliosphere. 



2.2. Coordinate systems 

We analyze ACE data in Geocentric Solar Ecliptic (GSE) sys- 
tem of reference {xgse, ycsE, icsE), where Xcse points from 
the Earth toward the Sun, ycsE is in the ecliptic plane and in the 
direction opposite to the planetary motion, and i.csE points to 
the north pole. However, Ulysses data are provided in the heli- 



M.S. Nakwacki et al.: Evolution of a MC from Sun to 5.4 AU 



3 



ographic Radial Tangential Normal (RTN) system of reference 
(R, T, N), in which R points from the Sun to the spacecraft, T 
is the cross product of the Sun's rotation unit vector (Cl) with R, 
and N completes the right-handed system (e.g., iFranz & Harped 
I2002h . 

In order to make an accurate comparison between the ob- 
servations of vector quantities made from both spacecraft and 
the orientation of the flux rope at both locations, we rotate all 
the vectorial data from Ulysses to the local GSE system of ACE 
(when this spacecraft is at the closest approach distance of the 
MC axis). We describe this transformation of coordinates in 
Appendix |A] This permits to compare magnetic field compo- 
nents in the same frame (Figs. I1I2I ). as well as to compare the 
orientation of the MC at both positions. 

We next define a local system of coordinat es linked to the 
cloud (i.e., the cloud frame. lEepping et aHl 19901) in order to bet- 
ter understand the cloud properties and to compare the results at 
both positions (such as the axial/azimuthal magnetic flux). The 
local axis direction of the MC defines Zdoud (with B;,cioud > 0). 
Since the speed of the cloud is mainly in the Sun-Earth direc- 
tion and is much larger than the spacecraft speed, which can be 
supposed to be at rest during the cloud observing time, we as- 
sume a rectilinear spacecraft trajectory in the cloud frame. The 
trajectory defines a direction d; so, we take ydoud in the direc- 
tion Zcioud X d and x^oud completes the right-handed orthonormal 
base (Xcioud, ydoud, Zcioud)- Thus, B.t.cioud, fi.v.cioud, fi;,cioud are the 
components of B in this new base. 

The cloud frame is especially useful when the impact pa- 
rameter, p (the minimum distance from the spacecraft to the 
cloud axis), is small compared to the MC radius (called R be- 
low). In particular, for p - and a MC described using a cylin- 
drical magnetic configuration, B{r) = B~{r)z + B^{r)^, we have 
Xcioud - f and ycioud - 4> after the spacecraft has crossed the MC 
axis. In this case, and for a cylindrical flux rope, the magnetic 
field data obtained by the spacecraft will show: Becloud = 0, a 
large and coherent variation of By.cioud (with a change of sign), 
and an intermediate and coherent variation of B;,cioud, from low 
values at one cloud edge, taking the largest value at its axis and 
returning to low values at the other edge (B;,cioud = is typically 
taken as the MC boundary). 

One possible procedure to estimate the flux rope orientation 
is the classical minimum variance (MV) method applied to the 
normalized series of magnetic field measurements within the es- 
timated boundaries of the MC dSonnerup & Cahilll[l967h . It was 
extensively used to estimate the orientation of MCs (see e.g., 
Lepping et al. 1990t Bothmer & SchwennI 19981 : iFarrugia et af] 
1999t iDasso et al J 120031: iGuUsano et al.ll2005h and it provides 
a good orientation estimation when p is small compared to R 
and if the i n/out bound m a gnetic fields are not significantly 
asymmetric. iGulisano et alJ (l2007h have tested the MV using 
a static cylindrical Lundquist's solution. They found a devia- 
tion of the axis orientation from the model of typically 3° for 
p being 30% of R. This deviation remains below 20° for p as 
high as 90% of R. Another method to find the MC orientation is 
called Simultaneous Fitting (SF). It minimizes a residual func- 
tion, which takes into account the distance between the observed 
time series of the magnetic field and a theoretical expression 
containing several free parameters, which include the angles for 
the flux rope orientation and some physical parameters associ- 
ated with the physical model assumed for the magnetic configu - 
ration in the cloud (e.g. iHidalgo et alJl2002l:lDasso et alJl20()3h . 



3. Modeling the magnetic cloud evolution 

3. 1 . Self-similar expansion 

The ev olution of a MC can be described with the model devel- 
oped bv lDemoulin et al.l(l2008l) . In this model, based on previous 
observations and theoretical considerations, a few basic hypoth- 
esis are introduced. Firstly, the MC dynamical evolution is split 
in two different motions: (i) a global one describing the position 
fcuit) - D{t)\cM of the center of mass (CM) with respect to a 
fixed heliospheric frame and (ii) an internal expansion where the 
elements of fluid are described with respect to the CM frame. 
Secondly, during the spacecraft crossing of the MC, the motion 
of the MC center is approximately a uniformly accelerated mo- 
tion and thus 



D{t) = Do + Voit - to) + a{t - tofll . 



(1) 



Thirdly, the cloud coordinate system, {Xdoud, f cloud, ^cUmd) de- 
fines the three principal directions of expansion. Fourthly, the 
expansion of the flux rope is self-similar with different expan- 
sion rates in each of the three cloud main axis. In the CM frame, 
this assumption implies that the position, r(f), of a element of 
fluid is described by 



r(f) = x(f) + y(t) ydoud + z{t) 

= Xq e(t) ±cloud + yo f(t) ydoud + ZQ g{t) Zdoud , 



(2) 
(3) 



where x{f),y{f),z{f) are the fluid coordinates from the CM refer- 
ence point at time f, and where xo,yQ,ZQ are the position coordi- 
nates taken at a reference time fo- The time functions e(f), f{t), 
and g{t), provide the specific time functions for the self-similar 
evolution. Finally, based on observat ions of diff'erent MCs at dif- 
ferent distances from the Sun (e.g., Liu et alJ 20051 : IWang et aP 
l2005at iLeitner et al1l2007t iGulisano et alJ l20IOh . we approxi- 
mate e(f) by the function: 



e{t) = (D{t)/Do)' , 



(4) 



and similar expressions for f{t) and g{t), simply replacing the 
exponent I by m and «, respectively, in order to permit an 
anisotropic expansion, and proton plasma beta (fip). 

From the conservation of mass we model the decay of the 
proton density as 



-(l+m+n) 



(5) 



From the kinematic self-similar expansion proposed before, and 
assuming an ideal evolution (i.e., non-dissipative, so that the 
magnetic flux across any material surface is conserved), the evo- 
lution of the magnetic components advected by the fluid is 



B 



.t.cloiid - B x,c\o\iAo^D I Da) 
By,cloud = By.cloudo(£>/£>0)"''^'" , (6) 
B;,cloiid = B~^c\oudo(D/Do) 

With the above hypothesis and neglecting the evolution of 
the spacecraft position during the MC observation, the observed 
velocity profile (Vy) along the direc tion \c^f of the center of 
mass velocity is expected to be (.Demouhn et al.ll2008l) : 



V, = -Vo - a{t - to) + Vo 
y2 

* -Vo + -f^^t - to) , 
Uo 



t-to 



Do/Vo + t-to 



(7) 
(8) 
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where y is the angle between Zdoud and Vcm and 



3.3. Global magnetic quantities using the Lundquist's model 



I sin^ y + n cos^ y . 



(9) 



For typical values of MCs we can linearize Eg. Q in r - fo and 
neglect the acceleration a (iPemoulin et al.l 120081 " This implies 
that the slope of the observed linear velocity profile provides 
information on the expansion rate of the flux rope in the two 
combined directions: Xdoud and idoud (since ( involves both / 
and n). 

3.2. Magnetic field 

Since MCs have low plasma /? (a state near to a force free field) 
and present flux rope signatures, its magnetic configuration is 
generally modeled using the cylindrical linear force-free field 

- Bq{J \{aQr)(f>+J(){aQr)i] (Lundquist 1950). If the expansion 
coefficients in the three main cloud axis {I, m, n) would be signifi- 
cantly different, then an initial Lundquist configuration would be 
strongly deformed. However, observations of the MC field con- 
figuration are approximately consistent with this magnetic con- 
figurat ion at different helio-distances, ranging from 0.3 to 5 AU 
(e.g., iBothmer & Schwennl[T99i: iLeitner et alJl2007h . So that 
we expect a small anisotropy on the expansion along different 
cloud directions (i.e., I x m ^ n). From observat ions of dif- 
ferent MCs at significantly different helio-distances dWang et alJ 
l2005bt ILeitner et alJ2007h and from observatio ns of the velocity 
profile slope from single s atellite observations (iDemoulin et alj 
l2008HGulisano et al.l20l"oh . it has been found that ( ^ 0.8. Since 
^ is a combination of / and n, which depends on y for each cloud, 
a systematic difference on / and « would be detected in a set of 
MCs with variable y angle. Such systematic variation of I and « 
with y was not found. 

The kinematic self-similar expansion given in the previous 
section combined with a non-dissipative regime, as expected for 
space plasmas, provide a prediction for the observed magnetic 
configuration during the transit of the MC. Then, assuming an 
initial Lundquist confi guration, the observations are modeled as 
dPemouUn et al.ll2008h : 



fi.v,cloud(0 
■Sv,cloud(0 
Bz,cloudit) 



Bo 



Jl[U{t)], 



Pit) f{t)git) 

Vc(f-fc)siny Bq 



Bq 



e(t)f{t) 



pit) eit)g{t) 
MUit)] , 



JdUit)]. 



(10) 
(11) 
(12) 



and proton plasma beta (B,,). where U(t) - ""P^'^i p(j'^ — 

y/iVcit - fc) sin y)^ 4- p^, Bq is the strength of the magnetic field 
and Q'o/2 is the twist of the magnetic field lines near the center, 
at time t = t^.By construction of the self-similar expansion, this 
magnetic field is divergence free at any time. 

Equations ( fT0lfT2l i have free parameters which are computed 
by fitting these equations to in situ observations, by minimizing 
a residual function and quantifying the square of the difference 
between the observed and the predicted value for the magnetic 
field components (i.e., a least square fit). This provides informa- 
tion on the observed flux rope, such as its orientation, its exten- 
sion and its magnetic flux. We call this method the Expansion 
Fitting (EE) method, and EFI method when isotropy {I - m - n) 
is assumed. 



Quantification of global magnetic quantities, such as the so- 
called ideal-MHD invariants, has been very useful to com- 
pare and associate MCs and their solar sources (se e, e.g., 
'M andrini et all 120051: lOasso et al.l l2005bt iDassol l2009l) . These 
MHD invariants are computed using a specific model of the MC 
magnetic configuration (Section lJSb . 

For the Lundquist's solution, the axial flux is 



z,Limd 



= 27r r 
Jo 



B- rdr - InJiiaoRo)- 



BnR, 



0«0 



O'O 



(13) 



where Ro is the flux rope radius at a reference time fo- The az- 
imuthal flux is 



^ ip,Lund 



Jo 



B^dr. 



(1 - Jo(aoRo)) 



BqLq 



(14) 



where Lq is the axial length of the flux rope at to. For a flux rope 
staying rooted to the Sun, L(t) is typically of the order of the 
distance to the Sun D(t). 

The relative mag netic helicity is (iDasso et al ] 120031: 
iNakwacki et~ai]|20081: 



Lund 



AnL 



A^B^rAr 



= 2;r U, 



■A- 



2JoJi \ LoBM: 



0-"^o 



ao^o / O'O 



(15) 



(16) 



The magnetic energy content is not an invariant in MCs. In 
order to compute its decay rate we simplify and assume that the 
MC expansion is isotropic with e{t) - f{t) - g(t). Then, the 
magnetic energy (Eiund) is computed as ( INakwacki et alj2008bh 



-'Lund 



(t) 



2nL r* 
2yuo Jo 



B^rdr 



e{t) 



aoRo 



^o^o^o 
Mo 



i(t()) 



2(f) 



(17) 



where fio is the magnetic permeability. Thus, for the isotropic 
expansion, the magnetic energy decays with time as e{t)^\ In 
the case of a small anisotropic expansion, a similar decay is ex- 
pected. 

3.4. Global magnetic quantities using the direct method 

We define below a method to estimate the global magnetic quan- 
tities directly from the observations. This method assumes firstly 
that the cross section is circular, secondly that there is symme- 
try of translation of B along the main axis of the flux rope, and 
finally that the impact parameter is low (so that. By ^ iB^ and 
r R! X). We also separate the time range covered by the cloud 
in two branches and proton plasma beta (J3p). (the in-bound/out- 
bound branches corresponding to the data before/a fter the clos- 
est app roach distance to the flux rope axis, see e.g. iDasso et"al1 
l2005al) . and we consider each branch in a separate way. We 
call this method DM-in/DM-out, depending on the branch (in- 
bound/out-bound) that is used. We define the accumulative mag- 
netic fluxes for the axial and azimuthal field components, in the 
in-bound branch as 



Fz.DM-in{x) - In I B,(x')x'dx' 

JXi„ 



y,DM~m 



(x) 



By(x') dx' , 



(18) 
(19) 
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Table 1. Timings (dd, hh:mm UT) of the interfaces for sub- 
structures inside the ICME, identified with numbered ticks in 
Figures [T]|4] 



Tick number 


Timing-ACE 


Timing-Ulysses 


Substructure 


1 


04, 11:00 


23, 13:30 


sheath 


2 


04, 14:30 


23, 22:00 


MC in-bound 


3 


05, 07:45 


26, 13:00 


MC out-bound 


4 


05, 20:30 


28, 00:00 


back 


5 


06, 02:30 


28, 09:00 





where Xin is the x value at the starting point of the flux rope 
and dx - dt. Next, from By(x) an d F^j)M-in(x) we o btain an 
expression for the magnetic helicity (iDasso et al.ll2006h 

HDM-in=2Ll F-,DM-m{x)By{x)xdx. (20) 

Finally, the magnetic energy is computed from the direct ob- 
servations of B{t) - \B(t)\, 

EoM-in^ — B^(x)xdx. (21) 

In an equivalent way, we define the same quantities for the 
out-bound branch. These equations are used to estimate these 
global quantities directly from in situ observations of the mag- 
netic field, using the components of the field and the integration 
variable (x) in the local MC frame. 

4. Results for the studied MC 

4. 1 . Common frame for data at ACE and Ulysses 

We first describe the in situ data of both spacecraft in GSE 
coordinates (defined at the time of ACE observations, see 
Appendix |A]i. Figure [1] shows the magnetic field and plasma 
ACE observations. The flux rope extends from ticks '2' to '4'. 
There, the magnetic field is stronger by a factor x 3 than in the 
surrounding S W and it has a coherent rotation. At ACE position 
the velocity profile is almost linear within the MC (Figure [T]i. 
The density is relatively important, up to x 20 cm"^, com- 
pared to the density present in the SW before the MC sheath, 
« 5 cm"^. In most of the MC the proton temperature is clearly 
lower than the expe cted temperature (red line ) in a mean SW 
with the same sp eed (iLopez & Freeman 1986: Elliott et al.l2005t 
lDemouhnll2009h : this is a classical property of MCs. 

Figure |2] shows Ulysses observations, in the same format as 
Figure [1] At Ulysses position, the linear velocity profile is still 
present before a strong shock on March 26, 22:20 UT (tick 'S' 
in Figure |2]i. The proton temperature is, as typically observed 
in MCs, well below the temperature expected for a typical SW 
at 5.4 AU from an extrapolation of the empirical law given by 
iLopez &" Freeman (1986). Conversely to ACE, the density is 
much lower in the MC than in the SW present before the MC 
sheath (* 0.1 - 0.2 cm"-' in the MC compared to » 0.5 - 1 cm""* 
in the SW). With mass conservation, this implies that the vol- 
ume expansion rate of the MC is much higher than the SW one. 
The magnetic field observed at Ulysses, has very significantly 
decreased with respect to the field at ACE (factor x 20), be- 
coming even weaker than the one present in the surrounding S W 



Table 2. Magnetic cloud orientation according to minimum vari- 
ance and simultaneous fitting. 



Model 


Parameter 


ACE 


Ulysses 


Min. Var 


e 


12° 


-29° 


Sim. Fit 


9 


-11° 


-14° 


Min. Var 


<P 


101° 


84° 


Sim. Fit 





115° 


81° 



(typically by a factor » 2). This is consistent with the important 
observed decrease in density. 

4.2. Magnetic cloud at ACE 

The definition of the MC boundaries is an important step in the 
analysis of a MC since the selected boundaries are affecting all 
the physical quantities related to the magnetic field. The MC 
boundaries are associated to discontinuities in the magnetic field 
because such discontinuities are formed in general at the bound- 
ary of two regions having different magnetic connectivities, such 
as the flux rope and its surrounding medium magnetically linked 
to the SW. 

The front boundary of MCs is typically better defined than 
the back (or rear) boundary. Such is the case in the analyzed 
MC at ACE, where there is a fast forward shock at tick '1', 
and a strong discontinuity of (pB (at tick '2' in Figure [TJ. 
Moreover, the magnetic field in front has a reverse (pB and, 
between ticks T and '2', B is fluctuating, a characteristic of 
MC sheaths. The density has also a discontinuity and is en- 
hanced just before tick '2', another characteristic of MC sheaths. 
Then, the MC front boundary is set at tick '2'. It is worth not- 
ing that the shock at tick ' 1 ' was previously identified in Wind 
spacecraft data at 11:05 UT as an ICME related shock (see 
[http://pwg.gsfc.nasa.gov/wind/currentJistIPS.htm). This inter- 
planetary fast forward shock -wave has been recent ly studied in 
a multispacecraft analysis bv lKoval & Szabol (l2010l) . For further 
details on the characteristics of shocks and their identification 
in the interplanetary med ium see Vinas & Scudder (1986) and 
iBerdichevskv et al.l(l2000h . 

The back boundary is also set at a discontinuity of the mag- 
netic field. As typical in MCs, there are several possibilities after 
18 UT on March 5 (see Figure[T]). However, (f>B and Ob have clear 
discontinuities, similar but weaker than the front one, at tick '4'. 
This is confirmed by a discontinuity in the density. The above 
boundaries are used to find the orientation of the flux rope both 
with the MV and the SF (Section IZSli. Indeed, the MC bound- 
aries are better defined in the MC frame since the axial and az- 
imuthal field components are separated (Figure [3). So we need 
to determine the MC frame. 

We first use the MV method to find the direction of the MC 
axis. The MV should be applied only to the flux rope, otherwise 
if part of the back is taken into account the directions given by 
the MV could be significantly bias (the back regi on is no longer 
part o f the flux rope at the observation time, iGulisano et al.l 
I2007h . Then, an iteration is needed starting with the first esti- 
mation of the flux-rope boundaries from the data, performing 
the MV analysis, then plotting the magnetic field in the cloud 
coordinates, and finally checking if the selected back boundary 
is correct with the accumulated azimuthal flux (Eq. (fT9]l). For the 
studied MC at ACE, the orientation found with the MV provides 
an almost vanishing accumulated flux at the back boundary se- 
lected above (Figure |3]l. Then, no iteration is needed, and the 
MV is providing a trustable orientation within the accuracy of 
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Fig. 1. In situ plasma and magnetic field of the ICME observed on March 1998 by ACE located at a; 1 AU from the Sun. From top 
to bottom: absolute value of the magnetic field {B - \B\, in nT), magnetic field vector orientation (GSE): latitude {9b) and longitude 
(0b), bulk velocity {V, in km s"') including in red the fitted straight line for the MC range (see Section l431 l. the expected (continuous 
red line) and observed (dots) proton temperature {Tp, in K), proton densit y (w,, , in cm""*), and proton plasma beta (J3p). Vertical lines 
mark different interfaces separating different plasma regions (see Section l4~2l for a description and Table[T]for timings). Horizontal 
dotted lines in 6b, <Pb, and fSp mark values at 0°, 180°, 1 as a reference, respectively. 



the method (Table|2]i. The small mean value of Be cloud indicates 
that the impact parameter is small (Figure |3]l. This implies that 
the orientation found by the MV method could diff'er from the 
real one by typically 10° dGuhsano et al.ll2b07) . 

The MC axis direction is also estimated with a standard si- 
multaneous fit (SF) of the Lundquist's solution to the observa- 
tions (see end of Section l272l i. The fitting minimizes the distance 
of the model to the observed magnetic field in a least square 
manner. As for the MV, it is important to take into account only 
the data in the range where the flux rope is present. The output of 
the SF method provides estimations for the MC frame vectors as 
the MV, and also the impact parameter and the physical param- 
eters (the free parameters of the Lundquist's solution). There is 
a significant difference in the MC axis orientations between the 
MV and the SF methods, larger in latitude (23°) than in longitude 
(14°). Moreover, with the SF results the accumulated azimuthal 
flux is significant at the back boundary (so the magnetic flux is 
not balanced), and there are no nearby significant discontinuities. 
We conclude that the MV method is providing a better approx- 



imation of the MC frame than the SF method in this particular 
MC. Since cp is close to 90°, it is a good assumption to consider 
that ACE crosses the flux rope front (or nose). Moreover since 6 
is small, the MC axis lies almost in the ecliptic plane. 

In summary, at ACE we select the MC boundaries as March 
4, 1998, at 14:30 UT and March 5, 1998, at 20:30. The front 
boundary differs by up to » 5h from previous studies. The larger 
difference is with the b ack boundary of the MC since it was set 
on M arch 06 at 03:1 5 dSkoug et al.l 12000) , 2:30 UT jLiu et al.l 
l2005h . and 06:30 UT (iDu et al.ll2007h . so significantlv later than 
our boundary set at tick '4'. Our boundaries take into account 
the extension of the flux rope when it was crossing ACE. In fact 
part of the MC characteristics are present after the back bound- 
ary '4' (e.g. strong and relatively coherent magnetic field, cooler 
temperature than expected), but some other quantities are closer 
to typical S W values such as Pp which is larger than 1 . This was 
previously found in other MCs ( Passo et al. 2006, 2007), and 
was called a back region. This type of region was interpreted as 
formed by a magnetic field and plasma belonging to the flux rope 
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Fig. 2. In situ plasma and magnetic field parameters of the ICME observed on March 1998 by Ulysses located at 5.4 AU from the 
Sun. The format is the same as for Figure[T] The magnetic field components are in the GSE frame defined at ACE (see Section lZ2] i. 
Vertical lines with the same reference number correspond to the same interfaces of substructures as at 1 AU. 



when it was close to the Sun, and later connected to the SW due 
to magnetic reconnection in the front of the flux rope, which in a 
low density plas ma as the SW c ould be more efficient due to the 
Hall effect (e.g.. [Morales et al.ll2005i) . So a back region has in- 
termediate properties between MC and SW since it is a mixture 
of the two. Indeed, this is the case at the time of ACE observa- 
tions since B;,cioud is fluctuating in the back region, while By^cioud 
retains more its coherency as in the flux rope (Figure O. 

4.3. Magnetic cloud at Ulysses 

At Ulysses the MC has a complex structure in most of the mea- 
sured parameters. The proton temperature is significantly below 
the expected temperature in the interval of time between tick 
T and about one day after tick '3' (Figure |2|. At the begin- 
ning of this time interval, the magnetic field has a coherent ro- 
tation, then later it is nearly constant up to a large discontinuity 
of the field strength (at tick 'S', ~ 0.4 days after tick '3'). Later, 
the magnetic field is stronger with a significant rotation. Then, 
the observations at Ulysses show a flux rope with characteris- 
tics significantly different from a "standard" flux rope such as 
observed at ACE. A magnetic field strength flatter than at ACE 



is expected if the flux rope transverse size increases much less 
rapidly than its length, since it implies that the axial component 
becomes dominant with increasing distance to the Sun so that the 
magnetic tension becomes r elatively weaker in the force balance 
(iDemoulin & Dassoll2009ah . The strong discontinuity followed 
by a strong field is more peculiar. Indeed, this MC was strongly 
overtaken by a faster structure, as sho wn by the velocity p anel of 
Figure |2] and previously identified bv lSkoug et al.l ( |2000|) . 

The proton temperature has a similar variation at both space- 
craft since it is becoming lower than the expected temperature 
approximately after the shock defined by the discontinuity of 
V. This discontinuity defines the position of tick T (Figure|2]i. 
The front boundary of the MC is less well defined than at ACE 
since there is no discontinuity. Still the behavior of the magnetic 
field is similar at both spacecraft as follows. Before tick '2', 6b 
is fluctuating while globally decreasing, while after tick '2' it 
is gradually increasing (with fluctuations) at both spacecraft. (pB 
has a global behavior similar to a step function at both space- 
craft, with nearly constant values both before tick ' 1 ' and after 
tick '2'. The main difference for (pB is a discontinuity at tick '2' 
at 1 AU, while (pB has a smooth transition at 5.4 AU. Then, we 
fix the beginning of the MC, tick '2', at the beginning of the pe- 



8 



M.S. Nakwacki et al.: Evolution of a MC from Sun to 5.4 AU 




> 



00 



> 



CO 



> 

N 



GO 




10 

5 

-5 
-10 



10 
5 


-5 
-10 



J L 








Un^jAiAiS^vW*^^^^ ^1 1 1 

\''^H\ 1 1 
~ . . . . 1 . 1 . . . 1 . . .1 1,1, 





04 



05 



06 



Time [day of March 1998] 



Fig. 3. Strength and components of the magnetic field vector at 1 AU in the cloud frame given by the MV method (see Section l4?2l l. 
Numbered ticks are the same as in Figure [T] The dashed blue lines are the Lundquist field model fitted with the EFI method 
(Section [32]i. The red thick line is the accumulative magnetic flux for the azimuthal component (in arbitrary units). 



riod where 6b starts increasing and (pB is nearly constant. This is 
confirmed in the MC frame, as B my ~ at tick '2' at both 1 
and 5.4 AU (Figs. [an. 

The back boundary is much more difficult to define since 
there are several possibilities, and the absence of a clear coher- 
ence between all the measured parameters. We tried several back 
boundaries, the MV and the SF methods, and use the iterative 
method described in Section 14.21 to check the selected bound- 
ary in the derived MC cloud frame. We found that the MC axis 
orientation can change by more than 20° and that the selected 
boundaries are not confirmed by the cancelation of the accumu- 
lated flux. 

The difficulties of both the MV and the SF methods are due 
to the large asymmetry of the MC when observed at Ulysses. 
Even normalizing the magnetic field strength is not sufficient 
since the out-bound branch is strongly distorted after the shock 
(tick 'S'). This is the consequence of the overtaking flow seen 
after tick '5' in Figure |2l The MC observati ons have character- 
istics comparable to the MHD simulations of lXiong et al.l ( l2007l 
l2009l) where they modeled the interaction of a flux rope with an- 



other faster one. At Ulysses, the internal shock has propagated 
nearly up to the MC center, so that most of the out bound branch 
is strongly distorted. This implies that both the MV and the SF 
methods cannot be used to find the MC orientation at Ulysses. 

In the MC frame deduced at ACE, the magnetic field com- 
ponents measured at Ulysses between ticks '2' and '3' have the 
expected behavior for a flux rope (Figure|4|, as follows. B.v,cioud 
is almost constant and small indicating a low impact parameter, 
fiv.cioud shows a clear rotation and B^.cioud is increasing from tick 
'2' to '3'. These are indications that the studied MC has not sig- 
nificantly changed its orientation (< 10°) between 1 and 5.4 AU, 
despite the presence of the SW overtaking the rear region of the 
flux rope. Therefore, for Ulysses we use below the MC frame 
found at ACE. 

The accumulated azimuthal flux is maximum after the inter- 
nal shock position (FigureUJi, an indication that the shock would 
have overtaken the flux rope center. However, we cannot trust 
the behavior of the accumulated azimuthal flux behind the shock 
since the orientation and the strength of the magnetic field are 
strongly modified. Indeed, all the magnetic field components are 
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Fig. 4. Strength and components of the magnetic field vector at 5.4 AU in the cloud frame given by the MV method applied at 1 AU 
on ACE data (see Section l431 l. The format is the same as for Figure[3] 



perturbed even before the shock, in the interval ['3','S']. Before 
tick '3', By,cioud is weak and it almost vanishes at tick '3', while 
fi,-,cioud is nearly maximum there, so we set the flux rope center 
at tick '3'. 

In summary, at Ulysses we initially select the ICME bound- 
aries, including the sheath, as from March 23 at 13:30 UT (tick 
'I' Uo March 28 at 09 :00 UT (tick '5'). The boundaries chosen 
bv iSkoug et aP (l2000l) were slightly different since they chose 
the range from March 24 at 02:00 UT to March 28 at 02:30 UT. 
The above detailed analysis of the magnetic field behavior in- 
dicates that the in-bound branch of the MC is from tick '2' to 
'3', with an out-bound branch from '3' to '4' strongly distorted. 
The in-bound extension is confirmed by a similar behavior of the 
velocity and the proton temperature at 1 and 5.4 AU. 



4.4. Impact parameter at ACE and Ulysses 

The agreement between the characteristics of the MC observed 
at ACE and Ulysses are clear indications that the same MC was 



observed. In this section we estimate the impact parameter when 
the MC is observed at each of the two spacecraft. 

Figure [3] shows that Bx.doud is small compared to B, as ex- 
pected when the impact parameter {p) is close to zero. From the 
m ean value of Brrinud, ^^ estimate pjR x 0.3 (see the method 
in lGulisano et al.ll2007l) . From the in-bound size Sacem =0.144 
AU and pjR - 0.3, we estimate that the radius Race ~ 0.15 AU, 
with a circular cross-section. Fitting the value of pIR, using the 
SF-EFI method, we also obtain pIR x 0.3. Another estimation 
of the impact parameter can be done using the approximation 
introduced in the simplest form of Eg. ( 31), < Bx^doud > I < 
B >a: \.2plR (iDemoulin & Dassol2009bh . assuming across sec- 
tion roughly circular We obtain p/R ^ 0.27 from this method. 
Thus, in conclusion, the impact parameter when ACE observes 
the cloud is low {p/R ~ 0.3) and very similar from its estimation 
using different proxies. 

From the mean value of B^.doud at Ulysses (see Figure|4]i, as 
done for ACE, we estimate pjR - 0.54. From S uiyj„ =0.55 AU 
and p/R - 0.54, we estimate that Rm^ x 0.65 AU. The above 
values of R at both spacecraft are in agreement with the values 
found from fitting the free parameters of the cylindrical expan- 
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sion Lundquist model (EFI method, see Section [3T2l Eqs. ( fTOl - 
fT2l) ). Again, using also the estimat ion of the impact parameter 
from the simplest approximation of lDemo ulin & Dasso (2009b) 
(< B e cloud > / < -S >~ 1.2p//?), we obtain from this method that 
p/R 0.18. Thus, in conclusion, the impact parameter when 
Ulysses observes the cloud is also low, with p/R in the range 
[0.2-0.5] using different proxies. 

Further arguments in favor of the association of the MC ob- 
served at ACE and Ulysses are given in the following two sec- 
tions with the timing and the agreement of the mean velocity, 
expansion rate, magnetic fluxes and helicity as deduced at the 
two locations. 



Table 3. Observed (ACE and Ulysses) mean values of proton 
density (n^) and magnetic field (B, B^ mv, and B, mv) for the 
in-bound branch of the magnetic cloud, for three different esti- 
mations of expansion rates (see Section |4~6] |. The predicted to 
observed ratio, at Ulysses, is shown in parenthesis. 



Quantity 


ACE 


Uly 


UlyPl 


UlyP2 


UlyP3 


1 






0.7 


0.7 


0.74 


m 






0.7 


1. 


1. 


np (cm"'') 


15.8 


0.12 


0.27 (2.2) 


0.17 (1.4) 


0.16(1.3) 


B(nT) 


10.3 


0.48 


0.88 (1.7) 


0.56(1.1) 


0.53 (1.0) 


ByMv (nT) 


4.1 


0.18 


0.51 (2.8) 


0.24 (1.3) 


0.22 (1.2) 


B-^Mv (nT) 


9.0 


0.42 


0.85 (2.0) 


0.51 (1.2) 


0.48 (1.2) 



4.5. Acceleration and Expansion 

The translation velocity of the flux rope at ACE is estimated, 
as typically done, as the mean value of the observed bulk 
speed during the observation range ['2','4']. It turns out to be 
Vace - -348 km/s. However, for Ulysses, because of the per- 
turbation on the out-bound branch, it is not possible to apply 
this classical procedure. Then, we compute a mean value of 
the observed speed only in a symmetric range near the center 
(tick '3'). We choose a range of 12 hours around '3', which 
gives Vuiysses = -351 km/s. So that at Ulysses, the MC trav- 
els slightly faster than at ACE, with a mean acceleration of 
<a>- IVuiysses - VAcsl/Af = 0.14 km/s/d, where Af = 21.2 days 
is the elapsed time between both centers. However, the small 
value obtained for the acceleration just indicates that it was al- 
most negligible during the transit from ACE to Ulysses and it 
implies a negligible contribution to the interpretation of the ob- 
served velocity profile as a proxy of the expansion of the flux 
rope (Eq. (|7]i, see Demoulin et al. 2008, for a justification). Then, 
we consider below that the MC is traveling from ACE to Ulysses 
with a constant velocity. 

We fit the velocity observations to the velocity model de- 
scribed in Sect. 13.11 (Eq. [8]). From the fitted curve (red line, 
panel 4 Fig. [1] and Fig.|2| we obtain the expansion coefficients 
(ace - 0.74 for ACE and ^y/v = 0.67 for Ulysses. These values 
are very similar at both heUo-distances and they are consistent 
with previous r esults found at 1 AU w ith the same method for a 
set of 26 MCs (iDemoulin et al.ll2008h . as well as the results ob- 
tained from a statistical analy sis of MCs or ICMEs observed at 
various distances of the Sun ( Bothmer & Schwemill 19941: IChenI 



ll996HLiu et ai1l2005HLeitner et alJl2007h " 



The angle, y, between the axis of the MC and the direction 
of motion defines the contribution of the axial and ortho-radial 
expansion rate to the observed V i(f), see Eq. (|9]). From the orien- 
tation given by MV at ACE, this angle is y = 101°. This implies 
that the expansion rate measured from v v is mainly due to the ex- 
pansion in a direction perpendicular to the cloud axis and, then, 
for the cloud analyzed here ( ^ I. 

The value of ^ observed at a given spacecraft, just provides 
the 'local' expansion rate, which corresponds to the expansion 
during the m situ observations. However, since (ace ~ (uiv, we 
assume that during the full travel between ACE and Ulysses the 
mean expansion rate occurred at a value between ( = 0.67 and 
( = 0.74. 

Assuming a self-similar expansion in the Sun-spacecraft di- 
rection (as in Eq. (|4|), we can link (without any assumptions 
on the cloud shape and symmetries for the flux rope) the size 
of the structure along the direction of the MC motion (vcm), at 
both helio-distances with the expected size at Ulysses given by 

SuLY,exp — S ACE,obs(5-4-Y ■ 



Because the out-bound branch of the MC at Ulysses is 
strongly distorted, we analyze the in-bound branch. Then, us- 
ing the mean velocity and the time duration between ticks 
'2' and '3', we find that Sacejh =0.144 AU. Then, from the 
observed value of (ace, we find an expected size at Ulysses 
of S uiy.exp-iii =0.498 AU (using the observed (my, we obtain 
^'uiyexp-in =0-444 AU). This implies an expected center time on 
March 26, 9:00 UT (and even an earlier time when using (uiy)- 
However, the cumulative azimuthal flux {Fy, Figure |4] panel 3) 
is smoothly and monotonically increasing beyond this time, with 
a strong discontinuity of By a bit later, on March 26, 13:00 UT. 
At the center of the flux rope we expect a glob al extreme of Fy 
(iDasso et al l l2006i l2007t iGulisano etani2010l) . and, in particu- 
lar, as it was observed when this same cloud was located at 1 AU 
(panel 3 of Figure[3]l. However, it is not the case in Ulysses at the 
expected time (9:00 UT), so that we decide to take the center at 
13:00 UT, the time when By started to be disturbed. This lack of 
exact agreement between the predicted and observed center po- 
sitions could be associated with the not exact alignment between 
ACE and Ulysses. 



4.6. Prediction of ttie mean plasma density and magnetic 
field at Ulysses 

The expected values of the proton density and magnetic field 
at Ulysses can be predicted using the observations at ACE and 
the expansion rates along the three directions (/, m, and «, see 
Eqs. (|5]|6])). The mean expansion rate, (, along the plasma flow 
can be estimated with a value between (ace and (uiy, which for 
the orientation of this cloud results to be / ^ (. The presence 
of bi-direct ional electrons sup ports the connectivity of this MC 
to the Sun (ISkoug et al.ll2000l) : then, the axial expansion rate is 
estimated as « a; 1, because the axial length needs to evolve as D 
in or der to keep the magnetic con nectivity of the MC to the Sun 
(e.g.. IDemoulin & Dassoll2009ah . For the third expansion rate, 
m, we have no direct observational constraint. 

Because after 'S' at Ulysses (see Figure |2]i the cloud is 
strongly perturbed, we compare ACE and Ulysses in the in- 
bound branch. Column 2 of Table |3] shows the mean values in- 
side the in-bound branch, for the proton density (np) and the 
magnetic field observed at ACE. Column 3 shows the same mean 
values, but now observed at Ulysses. Columns 4-6 show three 
predictions for these quantities at Ulysses. The field strength, 
B, is computed from By Mv and B~ mv, since B^.mv is depen- 
dent on the impact parameter which differs at both spacecraft 
(anyway |B,v,mvI << B so that including By ^v does not change 
significantly the results). All the predictions are made using 
n - I. Column 4 (UlyPl) shows the prediction at Ulysses us- 
ing I - {(ace + (uiv)/'2- = 0.7 and m = I (i.e., an isotropic ex- 
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Table 4. MHD quantities calculated according to each model 
(see Sect. 13.21 ). The first column indicates the model, the second 
shows the name of the global MHD quantities and their units, 
the next two columns show ACE and Ulysses results, and the 
last column shows the percentage of decay between ACE and 
Ulysses results. 



Model 


Parameter 


ACE 


Ulysses 


% of decay 


DM-in 


F;/MX 10^' 


1.2 


0.8 


33 


EFI 


F;/MX 10^' 


1.0 


0.9 


11 


DM-in 


Fy/M^ 10" 


2.7 


2.5 


7 


EFI 


F,,/Mx 10" 


2.2 


1.9 


14 


DM-in 


H/Mx^ lO"*^ 


-6.5 


-3.9 


40 


EFI 


H/Mx- 10-*2 


-2. 


-1.8 


10 


DM-in 


£■/ erg lO^** 


18 


4. 


78 
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pansion in the plane perpendicular to the cloud axis). Column 5 
(UlyP2) shows the prediction at Ulysses using I = 0.7 and m - 1, 
which corresponds to an expansion such that the cross section of 
the magnetic cloud is deformed toward an oblate shape in the 
plane perpendicular to the cloud axis, with the major axis per- 
pendic ular to the global flow speed (e.g., iDemoulin & Dassol 
l2009bl) . Column 6 (UlyP3) shows the prediction at Ulysses us- 
ing / = 0.74 and m = 1, which corresponds to an expansion 
in the direction of the plasma main flow as observed at ACE, 
emulating the case in which the expansion in this direction was 
similar to that observed at ACE almost all the time. 

The assumption I = 0.70 and m = / (UlyPl) predicts sig- 
nificantly larger values for all the quantities with respect to the 
observed ones (Table |3]l. However, the assumption / = 0.70 and 
m-\ (UlyP2) provides more realistic predictions. Furthermore, 
the assumption m = 1, combined with I - 0.1 A, gives predictions 
closer to the observations. Of course due to the lack of a perfect 
alignment we do not expect an exact matching even when the 
expansion is well modeled. 

Another possible approach is to compute /, m, n from the ob- 
served ratio (Ulysses with respect to ACE observations) of the 
mean values of iip, By^MV, and B- mv in the in-bound. We find 
I = 0.78, m - 1.04, and n - 1.06. The value of / is close to ^ace, 
measured independently from in situ velocity, and n is close to 
the expected value obtained with a flux-rope length proportional 
to the distance to the Sun. 

We conclude that the expansion of this cloud between ACE 
and Ulysses was such that / s; (ace and m ^ n ^ 1 . From this 
anisotropic expansion, and assuming a circular cross section for 
the cloud at ACE, we predict an oblate shape at Ulysses with an 
aspect ratio of the order of (5.4)'"^''^'^ « 1.55. If this anisotropic 
expansion is present also before the cloud reaches 1 AU, this 
aspect ratio could be even a bit larger. 

4. 7. Magnetic fluxes and helicity 

From fitting all the free parameters of the expanding model to 
the observations (method EFI), we compute the global magnetic 
quantities (Section 13.31 Eqs. ( I13H17| |). We also compute these 
quantities from the direct observations in the in-bound branch, 
using the direct method (see Section [341 Eqs. [TsU^TI) . We use 
a length L - 2 for ACE and L = 10.8 for Ulysses, because the 
MC is still connected to the Sun when observed at 1 AU and at 
5.4 AU. There is a good agreement between the magnetic fluxes 
and helicities found at ACE and at Ulysses, with a trend to find 
slightly lower values at Ulysses, i.e., a small decay of x 7-40 %, 



depending on the method to estimate them (Table HJi. We recall 
that the EFI method uses the full MC observations, so it is less 
accurate because of the strongly perturbed out-bound branch. 

The magnetic energy is not an MHD invariant. In fact, its 
decay, assuming a self-similar expansion with I - m = n, is 
predicted as e(f) ' (Eq. (fTTTi). For the MC studied here, we ex- 
pect a decay with a factor 5.4"'. This factor is in the range 
[0.19,0.31] for I in the range [0.7,1.]. In fact, the results of 
Section 14.61 indicate that the expansion is anisotropic [l x 0.7, 
m X n a: 1.). The computation of the energy decay with such 
anisotropic evolution would require a theoretical development 
which is outside the scope of this paper. Still, the energy decay 
is expected to be within the above range. Since the anisotropy 
in the coefficients /, m, n is relatively small, an approximation of 
the magnetic energy decay is obtained using a mean expansion 
of (/ + m + n)/3 = 0.9, which implies an energy decay x 0.22. 
From the last two rows of Table 31 the observed decay between 
1 AU and 5.4 AU is 0.22 and 0.26 for DM and EFI, respectively. 
There is an excellent agreement with the theoretically expected 
decay, even when we have simplified the analysis to I - m - n 
and cylindrical symmetry. 

5. Solar source of the MC 

5. 1 . Searching for tlie solar source 

The first step to determine the MC source on the Sun is to delimit 
the time at which the solar event could have happened. We com- 
pute the approximate transit time from Sun to Earth using the 
MC average velocity at ACE (see Sect. 14. 5b . Considering that the 
cloud has travelled 1 AU at a constant velocity of Vc ~ 350 km/s 
(where we neglect the acceleration, which is important only in 
the first stages of the CME ejection), we find t a; 1 AU/Vc ~ 5 
days. As the structure was observed by ACE starting on 4 March, 
we search for solar ejective events that occurred 5 days before, 
around 28 February + 1 day. 

From 28 February to 1 March, 1998, 5 numbered active re- 
gions (ARs) were present on the solar disk (see top panel of 
Fig. |5]l. Only very low X-ray class flares occurred in this pe- 
riod of time, most of them were class B (3 on 27 February, 2 on 
28 February, 6 on 1 March) and 3 reached class C on 1 March 
(see the X-ray light curve from the Geostationary Operational 
Environmental Satellites in http://www.solarmonitororg and the 
corresponding list of events). Two of the latter C-class events oc- 
curred in AR 8169 which was located at a S21W74 at the time 
of the flares. We have found no AR associated to the observed B 
class flares. 

We have also looked for the CMEs that occurred from 27 
February to 1 March in the catalogue of the Large Angle and 
Spectroscopic Coronagraph (LASCO, iBrueckner et al.l ri995l) 
on board the Solar and Heliospheric Observatory (SOHO, see 
[http://cdaw.gsfc.nasa.gov/CME/list). Most CMEs on those days 
had angular widths not larger than 70° and originated from the 
eastern solar limb, except for a halo and a partial halo CME. 

The halo CME first appeared in LASCO C2 on 27 February 
at 20:07 UT. This was a poor event inserted in the LASCO cata- 
logue after a revision on January 2006. The CME is clearly vis- 
ible only in LASCO C2 running difference images (see bottom 
panel of Fig.|5]l, in particular after a 22:00 UT when its front has 
already left the C2 field of view (LASCO C2 field of view from 
20:07 UT until 22:08 UT is only partial). The speed of the CME 
leading edge seen in LASCO/C2 and C3 is a 420 km/s from 
a linear fit, while a fit with a second order polynomia provides 
a 340 km/s at 20 Rq; these values are in agreement with the ve- 
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locity measured at the front of the MC at ACE (a; 385 km/s), tak- 
ing into account that the MC velocity is expected to be slightly 
modified by the interaction with the surrounding wind during its 
travel to 1 AU. Moreover, a velocity of « 340 km/s gives a travel 
time of ^ 4.5 days, so as expected. 

On the other hand, the partial halo CME first appears in 
LASCO C2 on 28 February at 12:48 UT. Its central position 
angle is 236° and its speed from a second order fitting is ^ 
225 km/s, which is too low considering the MC arrival time at 
ACE. Furthermore, from an analysis of images of the Extreme 
Ultraviolet Imaging Telescope (EIT Delaboudiniere et al. 1995) 
on board SOHO in 195 A the CME seems to be a backside 
event. Therefore, the halo CME is the candidate to be the MC 
solar counterpart. 

To find the source of the halo CME on the Sun, we ana- 
lyze FIT images obtained on 27 February starting 2 hours be- 
fore the halo CME appearance in LASCO C2. FIT was work- 
ing in CME watch mode at that time and, therefore, only im- 
ages in 195 A with half spatial resolution and with a temporal 
cadence of a; 15 min, or larger, are available. A sequence of 
4 images with full spatial resolution in all FIT spectral bands 
was taken at a 07:00 UT, 13:00 UT, 19:00 UT Furthermore, 
there is an extended data gap in FIT starting at around 20:00 
UT until around 22:00 UT. Considering that all events on 27 
February were of very low intensity (a B2.3 flare at ^ 18:20 
UT, a B1.5 flare at a 19:17 UT, and a B4.3 at a 22:55 UT af- 
ter the halo CME), the low temporal cadence and spatial res- 
olution of FIT images, and the existence of the data gap, we 
have not been able to unambiguously identify the CME source 
region. However, from these images and th ose from the Soft 
X-ray Telescope (SXT, iTsuneta et al.l 1199 lb . it is evident that 
2 of the five ARs present on the solar disk displayed signa- 
tures of activity (compare, in particular, the image at 20:17 UT 
with the previous and following one in SXT AlMg movie in 
|http://cdaw.gsfc.nasa.gov/CME7Hst/); these are AR 8171 and 
AR 8164 (see Fig. Is). 

From the two ARs that could be the source of the halo CME 
on 27 February at 20:07 UT, AR 8171 was located at S24E02 
which is an appropriate location for a solar region to be the 
source of a cloud observed at Earth. However, the magnetic flux 
in this AR is a 2.0 x 10^' Mx. This value is small when com- 
pared to the range for the axial MC flux measured at ACE (1. - 
1.2 X 10^' Mx) sin ce, in general, a MC axial flux is 10% of the 
AR magnetic flux (iLepping et al.|[T997l) . Furthermore, accord- 
ing to the distribution of the photospheric fi eld of the AR polari- 
ties (i.e., the shape of m agnetic tongues, see lLopez Fuentes et al.l 
l2000tlLuoni et alj201 ll and Fig.|5]) the magnetic field helicity in 
this AR is positive, which is opposite to the MC magnetic helic- 
ity sign. Finally, the leading polarity of AR 8171 is negative, and 
with a positive helicity this implies that the magnetic field com- 
ponent along the polarity inversion line (PIL) points from solar 
east to west. Since the magnetic field component along the PIL 
is related to the axial MC field component, this is not compatible 
with the MC axial field orientation at ACE that points to the solar 
east (FigurelTJ. Therefore, we conclude that AR 8171 cannot be 
the solar source region of the halo CME despite its appropriate 
location on the disk. 

The other possible CME source region is AR 8 1 64, located at 
N16W32 on 27 February at a 20:00 UT. This region is far from 
the central meridian, considering this location and a radial ejec- 
tion, one would expect that ACE would have crossed the ejected 
flux rope eastern leg; however, ACE crossed the MC front (see 
Section|4j2]i which implies that during the ejection the flux rope 



suffered a deflection towards the east; this probably occurred low 
in the corona as the CME is a halo. Concerning the AR magnetic 
flux, its value is a 10. x 10^' Mx, which is high enough to ex- 
plain the MC axial magnetic flux as we discussed previously. 

The magnetic helicity sign of AR 8164 is negative, as shown 
by the shape and evolution of its photospheric polarities in 
Figure |6] (the spatial organization of the magnetic tongues on 
23-25 February). This sign is also confirmed by the coronal field 
model of the region (see Sect l5.2b and agrees with the MC helic- 
ity sign. 

Conversely to AR 8171, the leading polarity in AR 8164 is 
positive, implying a magnetic field component pointing from 
west to east along the PIL. In this case, this direction is com- 
patible with the orientation of MC axial field. Furthermore, the 
PIL forms an angle of a 60° in the clockwise direction with the 
solar equator, while the MC axis lies almost on the ecliptic (see 
Sect. I4.2l i. This difference between the PIL on the Sun and the 
MC axial direction can be explained by a counter-clockwise ro- 
tation of the ejected flux rope, as expected, since its helicity is 
negative (iTorok & Kliemll2005HGreen et al.ll2007h . 

From the previous analysis, we conclude that AR 8 164 is the 
most plausible source of the halo CME on 27 February, 1998, 
which can be the counterpart of the MC observed at ACE on 4-5 
March. In the next section we compute the magnetic helicity of 
the AR before and after the ejection and its variation; this value 
is used as a proxy of the magnetic helicity carried away from the 
Sun by the CME. 

5.2. Physical properties of fhe solar source 

Using AR 8164 MDl magnetograms, we have extrapolated the 
observed photospheric line of sight component of the field to the 
corona under the linear (or constant a) force-free field assump- 
tion: V X B = aB. We have used a fast Fourier transform method 
as proposed by |Alissandrakis (1981 ) and the tr ansformation of 
coordinates discussed in lDemouhn et al.l (Il997h . The value of a 
is chosen so as to best fit the observed coronal loops at a given 
time. We need high spatial resolution images to identify indepen- 
dent loops, which are not available at times close enough before 
the ejection. Then, we have used the full spatial resolution im- 
ages obtained by EIT in 171 A on 27 February at 13:00 UT; 
this is the closest time to the event in which coronal loops are 
visible. The boundary conditions for the model are given by the 
MDl magnetogram at 12:48 UT on the same day. The value of 
a is deter mined through an i terative process that has been ex- 
plained in iGreen etal T (l2002l) . The value of a that best fits the 
observed loops is a = -9.4x10"^ Mm"' (Figure|7]i. 

Once the coronal model is determined, we c ompute the rel - 
ative coronal magnetic helicity, Hcon following iBergej ( Il985h . 
In particul a r, we u se a linearized version of the expression given 
bv iBerge^ (Il985 | . see his Eg. A23) a s ha s been done in pre - 
vious works bv iMandrini et all ([2005') and iLuoni et al.l (l2005h . 
Following this approach, the magnetic helicity content in the 
coronal field before the ejection is Hcor = -1 1.4 x 10^^ Mx^. 

When a flux rope is ejected from the Sun into the IP medium, 
it carries part of the magnetic helicity contained in the coronal 
field. Therefore, we have to compute the variation of the coro- 
nal magnetic helicity, by subtracting its value before and after 
an eruptive event, to compare this quantity to the corresponding 
one in the associated IP event. As done before, we search for an 
FIT image in 171 A with full spatial resolution after the ejec- 
tion in which loops could be visible. However, in this case, the 
AR is much closer to the east limb and projection effects, added 
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Fig. 5. (top) MDI full disk magnetogram on 28 February, 1998, 
at 00:00 UT (positive/negative magnetic field polarities are indi- 
cated in white/black colour). The 5 ARs present on the solar disk 
are indicated with their NOAA numbers. Notice the location of 
AR 8171 and AR 8164. The presence of magnetic tongues with 
a shape compatible with a positive magnetic helicity sign is evi- 
dent in AR 8171. (bottom) LASCO C2 running difference image 
showing the halo CME at 22:08 UT together with the closest in 
time EIT running difference image in 195 A. 




- • / 



23 Feb. 1998 - 12:48 UT 



25 Feb. 1998 - 12:48 UT' 




Fig. 6. MDI magnetograms of AR 8164 showing the photo- 
spheric evolution of its main polarities. The shape and evolution 
of the magnetic tongues indicates that the magnetic helicity of 
AR 8164 is negative. The magnetograms have been rotated to 
the central meridian position of the AR (positive/negative mag- 
netic field polarities are indicated in white/black color). The size 
of the field of view is the same in all panels. 



the AR magnetic field is decaying, its flux is lower than before 
the CME (x! 7.0 x 10^' Mx); therefore, is also lower, Hcor 
= -8.1 X 10"*^ Mx". The real value of the coronal magnetic he- 
licity after the CME should be even lower than the later one, as 
we expect that the field relaxes to a closer to potential state. To 
determine the range of variation for He,,,-, we also compute its 
value taking the lowest a value {a = -6.3 xlO"^ Mm ') that still 
gives a good fitting to the global shape of EIT brightness after 
the CME; in this case, Ha,r - -5.4 x 10"*^ Mx". Considering the 
two values determined for Hcor after the CME, we estimate that 
3.3 X lO'*^ Mx^ < \AH,or\ < 6.0 x lO'*^ Mxl 

5.3. Link with tine observed MC 

From estimations of the helicity content when the cloud was ob- 
served at ACE and Ulysses, using the EFI and DM methods (see 
Table H, we found 2 x lO'*^ j^^2 < |^^^| ^ iq42 j^^i^ 

which is fully consistent with the range found for the release of 
magnetic helicity in the corona during the CME eruption. 

A fraction of the total magnetic flux of AR 8164 (a; 10 x 
10^' Mx) is enough to account for the magnetic flux in the MC 
at ACE, F, lO^i Mx) + F, (a; 2 x 10"' Mx) a 3 x lO^' Mx. 

Thus, we have found qualitative and quantitative proofs that 
let us associate the halo CME observed by LASCO C2 on 27 
February, 1998, to its solar source region (AR 8 164) and to its in- 
terplanetary counterpart, the MC observed at ACE on 4-5 March, 
1998, and at Ulysses on 24-28 March, 1998. 



to the low intensity of the coronal structures, make it even more 
difficult to distinguish the shape of individual loops; as a result, 
a cannot be unambiguously determined, i.e. we can adjust the 
global shape of EIT brightness with more than one a value. We 
selected the EIT image at 01:20 UT on 28 February and, follow- 
ing a conservative approach, we determine a lower bound for 
the coronal magnetic helicity variation. We select the closest in 
time MDI magnetogram (at 01:36 UT on 28 February) and, us- 
ing the previously determined value for a, we compute Hear- As 



6. Summary and Conclusions 

We have studied a magnetic cloud which was observed in situ 
by two spacecraft (ACE and Ulysses) in an almost radial align- 
ment with the Sun 2° for latitude and a 6° for longitude) and 
significantly separated in distance (ACE at 1 AU and Ulysses 
at 5.4 AU). This is an uncommon geometrical scenario and it is 
very appropriate for multi-spacecraft analysis of MC evolution. 
In each of the spacecraft locations, we have analyzed the cloud 
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27 Feb. 1998 - EIT 171 13:00 UT - MDI 12:48 UT 




Fig. 7. (top) EIT image in 171 A at 13:00 UT on 27 February, 
1998, with two MDI isocontours overlaid (+100 G) shown with 
continuous lines (positive: white line, negative: black line). The 
MDI magnetogram was taken at 12:48 UT. (bottom) The same 
EIT image with more isocontours (+ 100, 500 G) shown with 
grey lines and computed field lines superimposed (black contin- 
uous lines). Both axes are measured in Mm in the local solar 
frame. 



in the local frame (attached to the flux-rope axis) and quantified 
magnetic fluxes, helicity, and energy, using an expansion model 
of an initial Lundquist field (EFI) and a direct method (DM), 
which permits the computation of globa l magnetic quantitie s di- 
rectly from the observed magnetic field jDasso et al.l l 2005al) . We 
also computed the local non-dimensional expansion rate (^) at 
ACE and at Ulysses from the ob served bulk velocity profiles (as 
defined in lDemouUn et al.ll2008h . 

We found close values of the normalized expansion rate 
along the solar radial direction, [,ace - 0.74 and ^f//,, - 0.67, as 
measured from the radial proton velocity at ACE and Ulysses, 
respectively. From the measured I^ace in the radial direction and 
assuming a self-similar expansion proportional to the solar dis- 



tance in the ortho-radial directions, we successfully predicted at 
Ulysses the values of the MC size, the mean values for density 
and magnetic field components from the values of these quanti- 
ties measured at ACE. 

Next, comparing observations at Ulysses with different mod- 
els for anisotropic expansions in the two directions that cannot 
be directly observed {m and n), we found that the expansion on 
the plane perpendicular to the cloud axis is larger than in the di- 
rection perpendicular to the radial direction from the Sun. Based 
on the quantification of this anisotropic expansion, we conclude 
that the initial isotropic structure at lAU will develop an oblate 
shape such that its aspect ratio would be a; 1.6 at 5.4 AU (i.e., 
the major axis ^ 60% larger than the minor one, with the major 
one perpendicular to the radial direction from the Sun). 

From a comparison of the transit time, axis orientation, mag- 
netic fluxes, and magnetic helicity, and considering all the solar 
sources inside a time window, we have also identified the pos- 
sible source at the Sun for this event, finding an agreement be- 
tween the amounts of magnetic fluxes and helicity, in consis- 
tence with a rough conservation of these so-called ideal-MHD 
invariants. 

In particular, we found that there is a small decay of the mag- 
netic fluxes and helicity between 1 and 5.4 AU, with a a; 10% 
of decay for F- and F, , and a decay of =s 10% for the magnetic 
helicity when the EFI method is used and a 40% when DM is 
used, respectively. These decays can be due to a possible ero- 
sion or pealing of the flux rope during its travel, for instance be- 
cause of magnetic reconnection with the surrounding SW (e.g., 
Dasso et al. '2006i). 

For a self-similar expansion and known expansion rates, it is 
possible to theoretically derive the decay of the magnetic energy 
during the travel of the flux rope in the SW. From the observed 
values of ^ and modeling the expansion rates in the other two 
directions (m and n), we predict its decay during the travel from 
1 AU to 5.4 AU. The measurements confirm this expected mag- 
netic energy decay (from (15 - 18) x 10^*^ erg to 4 x 10^*^ erg). 

Summarizing, in this work we validate for the first time that 
the local expansion rate (^) observed from the velocity profile 
can be used to make predictions of the decay of mass den- 
sity and magnetic quantities. From the comparison of detailed 
predictions and observations of the decays of these quantities, 
we provide empirical evidence about the quantification of the 
anisotropic expansion of magnetic clouds beyond Earth, to 5 
AUs. Finally, we quantify how much the so-called ideal-MHD 
invariants are conserved in flux ropes traveling in the solar wind. 
Then, this kind of combined studies, using multi-spacecraft tech- 
niques, is a powerful approach to improve our knowledge of the 
properties and evolution of magnetized plasma structures ejected 
from the Sun. 



Appendix A: Transformation of coordinates 
between ACE and Ulysses 

We describe below the transformation of coordinates between 
the natural system where the vector data of Ulysses are provided 
(R, T,N) to the GSE system at the time when the cloud was 
observed by ACE. Then, we provide a common frame to com- 
pare vector observations made at Ulysses and ACE (as shown in 
Figures |3} ID . 

A third co ordinate system, the H AE system (Heliocentric 
Aries Ecliptic, iFranz & Harpeij (l2002b ). is used since the loca- 
tion of both spacecraft are known in the HAE system. In this 
frame, Tjhae is normal to the ecliptic plane, and X^ae is positive 
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Yhae 



Fig. 8. Schema of the positions of each spacecraft at the observed times (left). Schema showing the relationship between the GSE 
and the HAE system of coordinates (right). 



towards the first point of Aries (from Earth to Sun at the vernal 
equinox, a; 21 March). At the moment in which the cloud was 
observed by ACE, it was located at a longitude ca - 163° and 
a latitude ta = 0° in the HAE system. When the MC was ob- 
served by Ulysses, this spacecraft was located at a longitude of 
eu = 157° and at a latitude of tu - 2° . 

The solar equator is defined as the plane normal to the Sun's 
rotation vector (SI) and it is inclined hy a x 7.25° from Zhae- 
In 2000, the solar equator plane intersected the ecliptic plane 
at a HAE longitude of ^ 75.6°. Then, the angle between the 
projection of £1 on the ecliptic is /? a; 14.4°. When we write R 
and n in the HAE system of coordinates, we obtain: 



R = cos(t£/) cos{eu) Xhae + cos(tj/) sin(ej/) Yhae 

+ sin(Tj/) Zha£ , (A.l) 

n = sin(a) cos(/?) Xhae - sin(a) sm(J3) Yhae 

+ cos(a)ZHAE- (A.2) 

To obtain T and N at Ulysses in HAE coordinates, we compute: 



n xR 

T = 

|QxR| 

N = Rxt. 



(A.3) 
(A.4) 



To go to the local GSE system of reference we do a last rotation 
around Zhae, which coincides with Zgse, in an angle 6 - 180°- 
eA (see Fig.|8]l. 
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